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ABSTRACT 

ASCA and ROSAT X-ray data are used to obtain two-dimensional maps of the gas temperature, 
pressure and specific entropy in the Triangulum Australis cluster of galaxies. We find that this hot 
(T e = 10.3 ±0.8 keV) system probably has a temperature peak (T e > 12 keV) at the cluster core, which 
approximately corresponds to the adiabatic relation. An underdense gas sector, found in the ROSAT 
cluster image eastward of the core, has a higher temperature than average at that radius. At this higher 
temperature, the gas pressure in this region is equal to that of the rest of the cluster at the same radius, 
but the specific entropy of this gas is significantly higher (although the temperature difference itself is only 
marginally significant). We speculate that the existence of this region of underdense high-entropy gas, 
as well as the adiabatic central temperature peak, indicate recent or ongoing heating of the intergalactic 
medium in this cluster. The most probable source of such heating is a subcluster merger, for which 
the hydrodynamic simulations predict a qualitatively similar temperature structure. We point out that 
entropy maps can provide a physically meaningful way of diagnosing merging clusters and comparing 
the predictions of merger simulations to the data. 

Subject headings: galaxies: clusters: individual (Triangulum Australis) — intergalactic medium — X- 
rays: galaxies 



1. INTRODUCTION 

Spatially resolved measurements of the gas temperature 
in clusters of galaxies can provide valuable information on 
the dynamical history of these systems, pointing to those 
clusters with recent or ongoing merger activity (see, for 
example, hydrodynamic merger simulations by Schindler 
& Miiller 1993; Roettiger, Burns, & Loken 1993; Evrard, 
Metzler, & Navarro 1996). Such measurements for hot 
clusters have become possible with the advent of ASCA 
with its combination of angular and spectral resolution 
and the wide energy band of 0.5-10 keV (Tanaka, Inoue, 
& Holt 1994). Two-dimensional gas temperature maps of 
several clusters obtained with ASCA using independent 
techniques were presented, e.g., in Arnaud et al. (1994), 
Markevitch et al. (1994), Markevitch (1996), Churazov et 
al. (1996), Honda et al. (1996), Markevitch & Vikhlinin 
(1996), and Henriksen & Markevitch (1996). The lat- 
ter authors detected a prominent asymmetric temperature 
pattern in the cluster A754, confirming with improved ac- 
curacy a similar ROSAT finding of Henry & Briel (1995). 
The temperature nonuniformity in A754, together with 
its complex X-ray image, indicated that it is undergoing a 
major merger. A temperature structure of this kind was 
also found in A3558, as well as in Coma. On the other 
hand, no significant asymmetric features were detected in 
other well-resolved clusters, A2256, A2319 and AWM7, 
except for the lower temperatures at the positions of the 
infalling subunits in the former two systems. A significant 
radial temperature decline at large radii was found with 
ASCA in several clusters; for some of them, most notably 
A2256, ROSAT PSPC provided an independent confirma- 



tion (refs. above). 

Here we report on the detection of a possible merger sig- 
nature in another cluster, Triangulum Australis. It is a rel- 
atively nearby (z = 0.051) bright, hot system, which was 
overlooked in the optical band due to its low Galactic lati- 
tude. It was first discovered as an X-ray source (McHardy 



et al. 1981). From an EXOSAT observation, 
art (1991) obtained an average temperature 
keV and an absorption column of 2.5 j 
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somewhat higher than the Galactic value of 1.25 x 10 
cm~ 2 (Heiles & Cleary 1979). The cluster has since been 
observed with ROSAT, which produced a high-quality im- 
age and an accurate measurement of the absorption col- 
umn, and with ASCA, which provided a gas tempera- 
ture map. These results are presented below. At the 
cluster's redshift, lh' 1 Mpc = 25'.4, using H = 100 h 
km s" 1 Mpc^ 1 . 



2. ROSAT RESULTS 

The cluster was observed with ROSAT PSPC for 7 ks 
in September 1992. We have used the method and code 
of Snowden et al. (1994) to obtain the flat-fielded images 
of the source in the energy bands R2-R7 (0.20-0.42-0.52- 
0.70-0.91-1.32-2.01 keV), free of non-cosmic background. 
The cosmic component of the background was then calcu- 
lated individually in each energy interval using the outer 
parts of these images, with the assumed accuracy of 5% 
(Markevitch & Vikhlinin 1996). Fluxes from these images 
(excluding the poorly calibrated band R3) were fitted to 
estimate the cluster temperature and the absorption col- 
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2 Quoted errors are 90% one-parameter intervals throughout the paper. 
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umn. Wc did not include estimates of the PSPC calibra- 
tion uncertainty to our analysis, since our primary inter- 
est in the ROSAT data is to obtain the image and the 
absorption column value, well-determined by PSPC. Our 
best-fit N H value is 1.42±°; 2 ^ x 10 21 cm- 2 ^ and T e = 7 
(4-17) keV, the large temperature error being due to the 
high Nh- The absorption column is consistent with the 
Galactic value and is marginally allowed by EXOSAT . 
We also attempted to constrain the spectrum in the in- 
ner r — 1'.5 circle (half the core radius) of the cluster, to 
see if there is any trace of a cooling flow to which PSPC is 
highly sensitive. We approximately corrected for the en- 
ergy dependence of the PSPC PSF (Hasinger et al. 1992), 
which for this small radius makes a noticeable difference 
in the spectral fit. Fixing Njj at its value for the whole 
cluster and using only 0.5-2 keV data, we obtain for this 
region T e = 8.2 (4.2-oo) keV. Freeing Njj and including 
lower energies, the temperature is T e =5.0 (2.4-oo) keV 
and Nh = 1.7 ± 0.5 x 10 21 cm~ 2 . Inclusion of a cooling 
flow component does not improve the fit significantly; the 
best-fit inflow rate is 7 (0-33) h~ 2 Mq yr _1 with a simi- 
lar best-fit absorption. The cluster image does not have a 
brightness peak characteristic of a strong cooling flow, and 
our ASCA spectra (below) do not require a cooling flow 
component either. We conclude that there is no indication 
of a significant cooling flow in the center. 

The cluster image in the 0.5-2 keV band is presented 
by contours in Fig. 1. Its X-ray brightness peak has co- 
ordinates a = 16 h 38 m 20 s , <5 = -64°21'22" (J2000). As- 
suming T e =10 keV and an iron abundance of 0.26 (aver- 
age ASCA values, see below) and the absorption column 
derived above, the 0.5-2 keV cluster luminosity within 
r = 1 h- 1 Mpc is 2.1 ± 0.1 x 10 44 ft," 2 erg s _1 . Fitting the 
brightness profile in the range of radii 0-25' by a /3-modcl 
of the form I x (r) cx (1 + r 2 /a 2 )~ 3/3+1 / 2 (Jones & Forman 
1984) and ignoring the asymmetry of the image, we ob- 
tained a x = 3'.5±0'.2 (0.143 hr Y Mpc) and = 0.63±0.02. 
We will use the derived Njj value and the ROSAT image in 
our ASCA analysis below. Note the apparent underdense 
sector east of the cluster center, which will be discussed 
below. At the off-center distances of 3-10', the surface 
brightness in this sector is about 2 times lower than that 
at the same radius on the opposite side of the cluster. 



3. ASCA ANALYSIS 

ASCA observed Triangulum Australis in March 1995 in 
two pointings with 11.3 ks and 6.7 ks CIS useful time, 
which were offset by 7' from one another. The SIS has 
low sensitivity at high energies relevant for this hot clus- 
ter, and was operated in 1-CCD mode which covers only 
ll'xll' area, so we limited most of our analysis to the GIS 
data. Using the energy band 1-11 keV, fixing the absorp- 
tion column at the ROSAT best- fit value and fitting both 
pointings simultaneously, we obtained an average cluster 
temperature of 10.3 ± 0.8 keV and an iron abundance of 
0.26 ± 0.10 (relative to Allen 1973). The EXOSAT tem- 
perature value (8.0jZi 3 keV) is somewhat lower; however, 
assuming their higher best-fit value of Nh, our tempera- 
ture becomes 8.9 keV, consistent with EXOSAT. 

For the reconstruction of the two-dimensional map of 
the cluster gas temperature, the method described in de- 



tail in Markevitch (1996, and references therein) has been 
used. It consists of simultaneous fitting of the temperature 
in all image regions, using a 0.5-2 keV ROSAT PSPC im- 
age as a brightness template and modeling the ASCA mir- 
ror scattering. We restricted our spatially resolved anal- 
ysis to the 2.5-11 keV energy band due to the poor PSF 
calibration at lower energies. (A temperature value for the 
whole cluster obtained from this energy band is 10. 8l^' 3 
keV.) The estimated systematic uncertainties of the back- 
ground normalization (20%, la), the PSF model's wings 
(15%) and core (5%) and the ASCA effective area (5%) 
were taken into account, as well as the uncertainty due to 
the relative displacement of the ASCA and ROSAT im- 
ages and the statistical error of the ROSAT image. The 
nearest bright X-ray source to Triangulum Australis is 3° 
away, so the possibility of stray light contamination is ex- 
cluded (Ishisaki 1996). 

The image was divided into 3 concentric annuli with 
r = — 3 — 12 — 23' centered on the brightness peak. The 
second and third annuli were divided into 5 and 4 sec- 
tors, respectively, with one of the sectors coinciding with 
the low-brightness region. The spectra from each region, 
each pointing and each detector were binned in 7-10 en- 
ergy intervals for adequate Gaussian statistics. A uniform 
temperature over each image region was assumed. We ob- 
tained Xmin = 149 for 224 d.o.f. for the simultaneous fit. 
The low reduced x 2 indicates that our conservative com- 
pound errors are, if anything, slightly overestimated. The 
confidence intervals were then estimated by Monte-Carlo 
simulations, including systematic errors. Figures 1 and 
2a show the resulting temperature map. Two sectors to 
the north and west of the outer ring which were poorly 
covered by ASCA pointings are not shown, but they were 
accounted for in the fit for completeness. There are two 
features in the map which will be discussed below — a sig- 
nificant temperature peak in the cluster core, and a higher 
temperature in the cluster low-brightness sector 2 com- 
pared to other regions at this radius. (The contributions 
to the total Xmin due to these regions are 24 and 14 for 
28 data bins, respectively.) The latter temperature dif- 
ference by itself is statistically insignificant; however, due 
to the coincidence with the depression of the gas density, 
the region's specific entropy (Fig. 2 c) is significantly dif- 
ferent from the average at this radius. The entropy values 
in sectors 2-6 could not be drawn from the same gaussian 
distribution at the 95% confidence, due to the single large 
deviation in region 2. 

We note that the ASCA data alone do not require as 
large a central temperature peak. As a consistency check, 
we performed a fit with the emission measure of the cen- 
tral region (relative to other regions) being an additional 
free parameter. Normally, the 0.5-2 keV relative bright- 
ness in different regions is predetermined by the ROSAT 
image in our technique. (Note that the ROSAT-ASCA 
relative normalization is not used.) The resulting best-fit 
relative brightness in the central annulus was 22% higher 
than that from ROSAT, and the central temperature was 
lower, with the temperatures in all other regions, as well 
as error bars, increasing (dotted symbols in Fig. la). How- 
ever, both quantities were consistent within their 90-95% 
confidence intervals with the temperature obtained above 
and the ROSAT brightness, respectively. Larger temper- 
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ature errors in such a fit are natural because the impor- 
tant piece of information provided by the ROSAT image is 
withdrawn. We have checked the possibility that ROSAT 
underestimates the relative emission measure in the r < 3' 
region by 20%, for example, due to a greater absorption 
toward the center (a factor of > 2 increase is needed). 
The ROSAT best-fit N H in the central r = 3' region is 
1.55lg gg x 10 21 cm~ 2 , far short of what is needed to pro- 
duce such effect. The PSPC finite angular resolution de- 
pletes the flux in the image peak only by ~ 2%, which is 
absorbed by the assumed ASCA PSF uncertainty. A cool- 
ing flow, which would make the ROSAT image inadequate 
for our purpose, is not required by ROSAT (above), nor 
does it make any improvement to either the GIS or the 
SIS fit. (The short SIS exposure covering the central re- 
gion does not provide useful constraints for the main hot 
component.) 

An optical image of Triangulum Australis shows that 
there is a cD galaxy 0'.3 from the X-ray brightness peak. 
If it has an active nucleus, it may add a nonthermal com- 
ponent to the central flux. However, a 4.85 GHz radio im- 
age (Condon et al. 1993) shows no source at the galaxy's 
position with an upper limit of about 30 mJy (two or- 
ders of magnitude below the luminosity of the cD galaxy 
NGC 1275 in the Perseus cluster). Therefore, a signifi- 
cant nonthermal flux from this galaxy is unlikely, unless 
it is a radio-quiet AGN, unusual in a cluster center. The 
best-fit normalization of an additional X-ray component 
with a typical power-law AGN spectrum is zero, although 
anything more exotic (for example, a heavily self-absorbed 
AGN with L x (2-10 keV)~ 2 x 10 43 h~ 2 erg s" 1 ) cannot 
be excluded, due to the ASCA's limited energy band. 

We also tried to fit the relative brightness in sector 2 as 
another free parameter, which practically did not change 
the temperature and the normalization in this region com- 
pared to the value fixed by ROSAT. Given all the above, 
we will assume that the central discrepancy is probably 
a statistical deviation, and our temperature measurement 
using the ROSAT image is correct, although it appears 
to be model-dependent. Our major conclusions will hold 
even without the central temperature peak. 

4. DISCUSSION 

We will now discuss the implications of our tempera- 
ture measurement. Interestingly, from the ROSAT im- 
age alone, which is rather circular except for the low- 
brightness sector, one might expect that this sector has 
a higher temperature, if the gas is in pressure equilibrium 
and the gravitational potential of the cluster is reason- 
ably smooth and symmetric. Pressure variations not sup- 
ported by gravity would disperse in roughly their sound 
crossing time (~3x 10 s h~ 1 yr for the low brightness sec- 
tor). Since this sector presumably contains lower density 
gas, it must have a higher temperature to maintain a pres- 
sure similar to the adjoining sectors. As our temperature 
map indicates, the gas in this region is indeed likely to 
be hotter. It is useful to compare the gas pressure and 
specific entropy in different regions, which is done in Fig- 
ure 2(6, c). Here, the entropy per particle is defined as 
As = s - s = §fcln[(T/T )(p/po) -2/3 ] ; the subscript 



denoting central values. The emission measure- weighted 
average density in each sector was crudely estimated as the 
value from the symmetric /3-model distribution multiplied 
by the ratio of the square root of the actual X-ray bright- 
ness to that in the model. Figure 2b indicates that the gas 
pressure is nearly constant at a given radius. Thus, the gas 
could be nearly in hydrostatic equilibrium with a symmet- 
ric gravitational potential. Note, however, that the gas in 
region 2 is less dense than the average in the middle annu- 
lus and would be expected to rise buoyantly on a timescale 
of ~ 10 9 hr 1 yr. Thus, the surface brightness distribution 
in Figure 1 is not consistent with a stable equilibrium dis- 
tribution, and is likely to have resulted from a transient 
hydrodynamic process. We note that the image can per- 
haps be explained as well by a very irregular underlying 
mass configuration, but that would certainly mean that 
the cluster is merging. 

A given gravitational potential may support different 
gas configurations, which are hydrostatic but at the same 
time differ by their histories and their entropy distribu- 
tions. The specific entropy As, shown in Fig. 2c for our 
data, is a useful diagnostic of shock heating or of any 
other nonadiabatic process. It may provide information on 
how the particular, observed gas distribution has emerged. 
For comparison, Figure 2 c also shows the entropy for an 
isothermal /3-model (arbitrarily choosing a temperature of 
8 keV oLthe three cooler "reference" sectors of the second 
annulusEI). This is the kind of distribution cluster gas is 
expected to approach if left to itself, particularly if thermal 
conduction is effective. Two gas regions significantly de- 
viate from this isothermal equilibrium model: the cluster 
core (if we accept that our central temperature is correct) 
and sector 2. Another notable detail is that the gas in the 
core has nearly the same specific entropy as the average 
in the surrounding annulus. One explanation of this adi- 
abatic temperature peak is that the core gas has recently 
mixed turbulently with the surrounding gas. The higher 
temperature in the core would then result from the adi- 
abatic compression required to bring the pressure in the 
core into hydrostatic equilibrium in a given gravitational 
potential. The second annulus and the outer sectors have 
similar temperatures and different entropies, which sug- 
gests that any such mixing was confined to the central 
region. This is not unreasonable since the core has only 
about 10% of the gas within the outer boundary of our 
second annulus. 

However, the entropy in region 2 (the low brightness 
sector) is higher than the average for that annulus or that 
of the cluster core. This could not be due to mixing; it 
suggests that the gas in this region underwent a recent 
heating event. For example, this might have been due 
to the passage of a shock from a subcluster merger (e.g., 
Schindler & Muller 1993; Pearce, Thomas, & Couchman 
1994), or this gas may have been expelled from the core 
where the heating occurred. Presumably, this gas then 
expanded adiabatically into pressure equilibrium with the 
surrounding gas. Note that because the pressure has al- 
ready equalized but the entropy in this region is still high, 
the event should have happened about 10 9 h^ 1 yr ago. 
The same event may have heated the central gas suffi- 
ciently to produce a radially decreasing entropy gradient, 
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after which convective mixing would act to eliminate this 
gradient, producing a distribution with similar specific en- 
tropy in the core and the surrounding gas. One can es- 
timate the amount of energy which should have been in- 
jected into the gas residing in the core and region 2, as- 
suming that the cluster was initially isothermal at 8 keV, 
as Q w TAS - (3 - 6) x 10 61 h^l 2 erg (approximately 
equally distributed between the core and sector 2). To 
estimate the heating which might occur in a merger, we 
consider the kinetic energy acquired by the gas in two sub- 
clusters each with a total mass of 6 x 10 14 M Q inf ailing 
from a distance of 5 to 1 Mpc; this energy is about 10 times 
greater than required. Therefore, a merger of subclus- 
ters is a reasonable source of such heating. Hydrodynamic 
simulations predict that a head-on merger of subclusters 
should produce a hot core and an asymmetric temperature 
structure not unlike the observed one. 

Other sources of heating might include supernova-driven 
galactic winds (e.g., Evrard, Metzler, & Navarro 1996). 
However, the present day supernova rate expected for the 
stellar population in the core of this cluster is much too 
small to produce enough heating on the timescale limited 
by the possible lifetime of the observed temperature struc- 
ture. 

5. SUMMARY 



Using ASCA and ROSAT data, we have reconstructed 
two-dimensional distributions of the gas temperature, 
pressure and entropy in the Triangulum Australis clus- 
ter. We detect a significant temperature peak in the 
cluster core (although its existence is somewhat model- 
dependent), and a likely temperature increase in the sector 
coincident with the underdense region of the gas found in 
the cluster X-ray image. The underdense region is in pres- 
sure equilibrium with the gas at the same radius but has a 
higher specific entropy. We interpret the obtained entropy 
distribution as evidence of the ongoing or recent heating 
of the intracluster medium, perhaps by shocks from a sub- 
cluster merger. A temperature map with higher angular 
resolution to look into the structure of the cluster core, 
as well as covering the outermost cluster regions, would 
be useful for understanding the processes responsible for 
heating of the gas in this object, and perhaps in clusters 
in general. The entropy distribution is a useful diagnos- 
tic of shock heating or of any other nonadiabatic process. 
Therefore, it would be very helpful if simulations of cluster 
hydrodynamics were analyzed to produce the gas pressure 
and entropy maps, to enable a quantitative and physically 
meaningful comparison to the data. 

We thank the referee, A. C. Edge, for useful criticism. 
We gratefully acknowledge the support from NASA grants 
NAG5-2526 and NAG5-1891. 
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Fig. 1. — Temperature map of Triangulum Australis. Con- 
tours show the ROSAT PSPC surface brightness in the 0.5-2 
keV band (smoothed with a Gaussian of variable width for clar- 
ity). The contours are logarithmically spaced with a factor of 
s/2. The grayscale shows ASCA temperatures (darker is hot- 
ter). Regions (whole inner circle and two outer annuli divided 
into 5 and 4 sectors, respectively) are numbered and their tem- 
peratures are shown in Fig. 2a. The two sectors of the outer 
annulus which were poorly covered by ASCA are not shown. 
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Fig. 2. — Temperature (a), gas pressure (b) and specific en- 
tropy relative to the central value (c) in the cluster regions 
numbered in Fig. 1. Errors are 90%. The dotted symbols in 
panel (a) correspond to a fit with the relative central brightness 
being a free parameter (see the text). The solid horizontal lines 
correspond to the average values over each annulus, excluding 
sectors 2 and 3 for the second annulus. The dashed line in 
panel (c) corresponds to a symmetric isothermal /3-model with 
T e = 8 keV (the average temperature in our "reference" sec- 
tors 4, 5 and 6). The values of As used for an estimate of the 
amount of energy injected into the gas (see text) are shown by 
shaded bars. 



